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Abstract 

We show that the recent WMAP data can impose a generous upper bound 
on the effective masses of electron, muon and tau neutrinos defined in the 
kinematic measurements: (m)^ + (m)'f l + (m) 2 . = m\ +m| + m| < 0.5 eV 2 , or 
(m) a < 0.71 eV (for a = e,fj,,r). When current neutrino oscillation data are 
taken into account, we obtain (m) e < 0.24 eV and (m)^ ~ {m) T < 0.24 eV. 
Thus there is no hope to kinematically detect (m)^ and (m) T in any realistic 
experiments. 

PACS number(s): 14.60.Pq, 13.10.+q, 25.30.Pt 
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Thanks to the recent Super-Kamiokande [1], K2K [2], SNO [3] and KamLAND [4] exper- 
iments, we are now convinced that the deficit of atmospheric neutrinos and the deficit of 
solar v e neutrinos are both due to neutrino oscillations, a quantum phenomenon which can 
naturally happen if neutrinos are massive and lepton flavors are mixed. Furthermore, current 
experimental data indicate that solar and atmospheric neutrino oscillations are dominated 
respectively by v e — > and — > v T transitions. The neutrino mass-squared differences 
associated with solar and atmospheric neutrino oscillations are thus defined as 



where rrii (for % = 1,2,3) denote the physical masses of three neutrinos. Although a strong 
hierarchy between Am^ and Am^ tm has been observed, the absolute values of mi, m 2 and 
m 3 remain unknown. The kinematic limits on the effective masses of electron, muon and tau 
neutrinos can be obtained from the tritium /3-decay — ^He + e~ + V e , the tt + — > ji + + 
decay and the r — > 5n + v T (or r — > 3tt + v T ) decay, respectively. Today's results are [5] 



One can see that the experimental sensitivity for (m)^ is more than four orders of magni- 
tude smaller than that for (m) e , and the experimental sensitivity for (m) T is two orders of 
magnitude lower than that for (m) M [6]. Is there any hope to detect (m)^ and (m) T or to 
constrain them to a meaningful level of sensitivity? The answer to this question relies on 
how small (m)^ and (m) T are. 

The main purpose of this short note is to calculate (m) e , (m)^ and (m) T with the help of 
the recent WMAP data [7] and neutrino oscillation data. While (m) e has been extensively 
studied in the literature [8], a detailed analysis of (m) M and (m) T has not been done. It 
is therefore worthwhile to work out the upper bounds on (m) M and (m) T from currently 
available data, such that one may definitely answer the afore-raised question. Independent 
of the neutrino oscillation data, a generous upper limit (m) 2 e + (w) 2 -\- (m)^ = m\+m\+m\ < 
0.5 eV 2 can be achieved from the WMAP observation. Hence (m) a < 0.71 eV holds for 
a = e,/i or r. Such an upper bound for (m) a will be reduced by a factor of three, if 
current neutrino oscillation data are taken into account. In particular, m (m) T is a 

good approximation and both of them are insensitive to the Dirac CP- violating phase of 
the lepton flavor mixing matrix. We conclude that there is no hope to kinematically detect 
the effective masses of muon and tau neutrinos. 

Direct neutrino mass measurements are based on the analysis of the kinematics of charged 
particles produced together with neutrino flavor eigenstates \u a ) (for a = e, //, r), which are 
superpositions of neutrino mass eigenstates \vi) (for % = 1,2,3): 




(1) 



(m) e < 3 eV , 
(m) M < 0.19 MeV , 
(m) T < 18.2 MeV . 



(2) 




(3) 
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The unitary matrix V is just the lepton flavor mixing matrix. The effective masses of 
electron, muon and tau neutrinos in the kinematic measurements can then be defined [6]: 

(m)l = \V al \ 2 m\ + \V a2 \ 2 m 2 2 + \V a3 \ 2 m 2 , (4) 

for a = e, \i and r. The unitarity of V leads straightforwardly to a simple sum rule between 
(m) 2 and m 2 : 

(m)l + (m) 2 + (m)l = m\ + m\ + m\ . (5) 

Note that this sum rule allows us to derive an upper bound on (m) 2 independent of any 
neutrino oscillation data. This point can clearly be seen from the inequality 1 

m\ + ml + m 2 < (m ± + m 2 + m 3 ) 2 , (6) 

in which the sum of three neutrino masses has well be constrained by the recent WMAP 
data [7]: m\ + m 2 + m 3 < 0.71 eV at the 95% confidence level. Therefore, 

(m) 2 e + (m)l + (m) 2 T < 0.50 eV 2 . (7) 

This generous upper bound implies that (m) 2 a < 0.50 eV 2 or (m) a < 0.71 eV holds for 
a = e, fi and r. Two comments are then in order. 

1. The cosmological upper bound of (m)^ is more than five orders of magnitude smaller 
than its kinematic upper bound given in Eq. (2). In comparison, the upper limit of 
(m) T set by the WMAP data is more than seven orders of magnitude smaller than its 
kinematic upper limit. It seems hopeless to improve the sensitivity of the kinematic 
measurements to the level of 0.71 eV. 

2. The cosmological upper bound of (m) e is about four times smaller than its kinematic 
upper bound given in Eq. (2). The former may be accessible in the future KATRIN 
experiment [11], whose sensitivity is expected to be about 0.35 eV. 

If current data on neutrino oscillations are taken into account, however, more stringent 
upper limits can be obtained for the effective neutrino masses (m) e , and (m) T . 

To see how (m) a may be related to the parameters of neutrino oscillations, we make use 
of Eq. (1) to express mi and m 2 in terms of m 3 , Am 2 un and Am 2 tm . The results are 

m 1 = yfm%+ pAm 2 atm + qAm 2 nn , 

m 2 = \]ml+pAml tm , (8) 

where p = ±1 and q = ±1 stand for four possible patterns of the neutrino mass spectrum. 
For example, p = q = +1 corresponds to m\ > m 2 > m 3 , and so on. The present solar 



x It is worthwhile to mention that the sum m 2 = m\ + m 2 , + m 2 is a crucial quantity in the thermal 
leptogenesis mechanism [9], because it controls an important class of washout processes [10]. 
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neutrino oscillation data favor q = — 1 or mi < vn 2 [12], but the sign of p has not been fixed. 
Substituting Eq. (8) into Eq. (4), we obtain 



(m) a = ^/ml + p(l-\V a3 \ 2 )Aml tm + q\V al \2Aml un . (9) 

Taking account of IK3I 2 1 [13] and Am^ <C Am 2 tm , we arrive at the following approxi- 
mation for (m) e : 

(m) e « y/ml + pAmL . (10) 

Taking account of the observed (nearly) maximal mixing factor of atmospheric neutrino 
oscillations [1], which is equivalent to |V^ 3 | ~ |T^- 3 | ~ l/y/2 for \V e3 \ <C 1, we obtain 



« ym 2 + |Am 2 tm , 

Hr ~ ^3 + | Am atm ■ (H) 

One can see that (m)^ ~ (m) T is a natural consequence of current neutrino oscillation data. 
In addition, Eqs. (10) and (11) tell us that (m) e is slightly larger than (m)^ and (m) T for 
p = +1 or m 2 > m 3 ; and it is slightly smaller than (m) M and (m) T for p = — 1 or m 2 < m 3 . 

In general, one may analyze the correlation between (m) a and m 3 by use of Eq. (9), 
only if \V a i\ and \V a3 \ are already measured. Since the mixing angles of solar, atmospheric 
and CHOOZ (or Palo Verde) reactor [13] neutrino oscillations read as 

sin 2 2£ sun = 4|\/ el | 2 |\/ e2 | 2 , 
sin 2 2^ atm = 4|V^ 3 | 2 (l - |V^ 3 | 2 ) , 

sin 2 2# chz = 4|V e3 | 2 (l-|V; 3 | 2 ) , (12) 



we reversely obtain [14] 



V el I = ^=\/ cos 2 6 chz + V cos 4 6 chz - sin 2 26 sun , 



1 
1 

\V e3 \ = sin9 chz , 
|V^ 3 | = sin^atm , 



V e2 \ = -^=\/cos 2 9 chz - v/ cos 4 9 chz - sin 2 20 sun , 



\V t3 \ = \j cos 2 6 chz - sin 2 atm . (13) 

The other four matrix elements of V (i.e., |V^i|, |V^ 2 |, and IV^I) are unable to be 

determined from the afore-mentioned neutrino oscillation experiments. They can be derived 
from Eq. (13), however, if the Dirac CP- violating phase in the standard parametrization 
of V [15] is taken into account. The explicit expressions of |V^i|, (V^l, \V T i\ and \V T 2\ are 
given in Appendix A. We see that (m)^ and (m) T depend on the Dirac phase 5 in the chosen 
parametrization of V. However, the sensitivity of (m)^ or (m) T to 5 is negligibly weak. The 
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reason is simply that the contribution of cos 5 to |V^i| (or |V^i|) is suppressed by |14 3 | or 
sin^ c h z , and the contribution of |V^i| (or \V T i\) to (m)^ (or (m) T ) is further suppressed by 
Amg Un . It is therefore hopeless to probe the CP-violating phase 5, even if (m)^ and (m) T 
could be measured. 

To numerically determine the upper bound of (m) a , we need first of all work out the 
upper limit of m 3 set by the WMAP and neutrino oscillation data. In FIG. 1, we show the 
dependence of mi + vri2 + 777 3 on 777,3, where the best-fit values Am^ = 7.3 x 10~ 5 eV 2 and 
Am 2 tm = 2.5 x 1CT 3 eV 2 [12] have typically been input. Note that only the m 1 < m 2 case, 
which is supported by current solar neutrino oscillation data, is taken into account. For the 
m 2 < m 3 case, m 3 has an lower bound m 3 > \J Am 2 tm + Am 2 un ~ 0.051 eV; but for the 
m 2 > m 3 case, even m 3 = is allowed (inverted hierarchy). We see that these two cases 
become indistinguishable for m 3 > 0.2 eV, implying the near degeneracy of three neutrino 
masses. Once the WMAP limit mi + m 2 + m 3 < 0.71 eV is included, we immediately get 
m 3 < 0.24 eV. Then we have m 8 < 0.24 eV for i — 1,2 or 3. This result is apparently 
consistent with those obtained in Ref. [16]. 

Next we evaluate (m) Q by using the best-fit values sun 33° and 9 atm 45° [12] in 
addition to taking 6* c h z ~ 5° as a typical input, which is compatible with sin 2 29 c \ lz < 0.1 
extracted from the CHOOZ reactor neutrino experiment [13] 2 . Our numerical results for 
(m) e , (m) M and (m) r changing with m 3 are shown in FIG. 2. Once again, it is impossible 
to distinguish between the m 2 < m 3 case and the m 2 > m 3 case for m 3 > 0.2 eV. We find 
that the dependence of (m) a on m 3 is very similar to the dependence of mi + m 2 + m 3 on 
777,3. As expected in Eqs. (10) and (11), (m) e > (m) M (m) T holds for m 3 < m 2 (curves a 
and b in FIG. 2); and (m) e < (m) M ~ (m) T holds for m 3 > m 2 (curves c and d in FIG. 2). 
In view of the upper limit m 3 < 0.24 eV obtained above, we arrive at (m) a < 0.24 eV for 
a = e, ix or r. This upper bound is suppressed by a factor of three, compared to the upper 
bound obtained from Eq. (7) which is independent of the neutrino oscillation data. 

The result (m)^ ~ (m) T < 0.24 eV implies that there is no hope to kinematically detect 
the effective masses of muon and tau neutrinos. As the WMAP upper bound is in general 
valid for a sum of the masses of all possible relativistic neutrinos (no matter whether they are 
active or sterile), it seems unlikely to loosen the upper limit of (m) a obtained above in the 
assumption of only active neutrinos. Therefore, the kinematic measurements of (m)^ and 
(m) T have little chance to reveal the existence of any exotic neutral particles with masses 
much larger than the light neutrino masses. 

The author is grateful to W.L. Guo for some discussions. This work was supported in 
part by National Natural Science Foundation of China. 



2 Note that # c h z ~ 5° is also favored in a number of phenomenological models of lepton mass 
matrices. See Ref. [17] for a review with extensive references. 
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APPENDIX A 



This 



ent 



Appendix is devoted to the calculation of \V^i\, |V^ 2 |, and IK2I, which are 
unrestricted by current neutrino oscillation data. Taking account of the Dirac CP- 
c phase in the standard parametrization of V [151, we find 
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of |V^i], I V^ 2 1 , IK1I and | VV 2 | in terms of 9 sun , atm , chz and 5 can 



The explicit expressions of |V^i], (V^l, a 
then be obtained from Eq. (13) and Eq. (Al). 

cos 2 chz - sin 2 atm ( 1 + sin 2 6 chz ) , 

— ^ '- Vcos 4 chz - sin 2 20 sun 

2 cos 4 chz 



cos 2 atm 



+ 
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1/2 



(A2) 
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These results are useful for a numerical analysis of V by use 
on neutrino oscillations, if 5 is allowed to vary from to n. 



(A5) 

of current experimental data 



6 



REFERENCES 



[1] Y. Fukuda et ai, Phys. Rev. Lett. 81, 1562 (1998); ibid. 81, 4279 (1998); http://www- 
sk.icrr.u-tokyo.ac.jp/dpc/sk/. For a review, see: C.K. Jung, C. McGrew, T. Kajita, and 
T. Mann, Ann. Rev. Nucl. Part. Sci. 51, 451 (2001). 

[2] K2K Collaboration, M.H. Ahn et ai, Phys. Rev. Lett. 90, 041801 (2003). 

[3] SNO Collaboration, Q.R. Ahmad et ai, Phys. Rev. Lett. 89, 011301 (2002); Phys. Rev. 
Lett. 89, 011302 (2002). 

[4] KamLAND Collaboration, K. Eguchi et ai, Phys. Rev. Lett. 90, 021802 (2003). 

[5] Particle Data Group, K. Hagiwara et ai, Phys. Rev. D 66, 010001 (2002). 

[6] P. Vogel and A. Piepke, in Ref. [5]. 

[7] C.L. Bennett et ai, astro-ph/0302207; D.N. Spergel et ai, astro-ph/0302209. 
[8] For a recent review with extensive references, see: S.M. Bilenky, C. Giunti, J. A. Gri- 
fols, and E. Masso, hep-ph/0211462. See, also, Y. Farzan and A.Yu. Smirnov, hep- 
ph/0211341 (to appear in Phys. Lett. B). 
[9] M. Fukugita and T. Yanagida, Phys. Lett. B 174, 45 (1986). 
[10] W. Buchmiiller, P. Di Bari, and M. Plumacher, hep-ph/0302092. 
[11] KATRIN Collaboration, A. Osipowicz et al, hep-ex/0109033. 

[12] V. Barger and D. Marfatia, hep-ph/0212126; G.L. Fogli et ai, hep-ph/0212127; M. 
Maltoni, T. Schwetz, and J.W.F. Valle, hep-ph/0212129; J.N. Bahcall, M.C. Gonzalez- 
Garcia, and C. Pena-Garay, hep-ph/0212147; P. Aliani, V. Antonelli, M. Picariello, 
and E. Torrente-Lujan, hep-ph/0212212; P.C. de Holanda and A.Yu. Smirnov, hep- 
ph/0212270; G.L. Fogli, E. Lisi, A. Marrone, and D. Montanino, hep-ph/0303064. 

[13] CHOOZ Collaboration, M. Apollonio et ai, Phys. Lett. B 420, 397 (1998); Palo Verde 
Collaboration, F. Boehm et ai, Phys. Rev. Lett. 84, 3764 (2000). 

[14] Z.Z. Xing, Phys. Rev. D 65, 077302 (2002); W.L. Guo and Z.Z. Xing, Phys. Rev. D 67, 
053002 (2003). 

[15] See, e.g., H. Fritzsch and Z.Z. Xing, Phys. Lett. B 517, 363 (2001); Z.Z. Xing, Phys. 

Rev. D 65, 113010 (2002). 
[16] A. Pierce and H. Murayama, hep-ph/0302131; C. Giunti, hep-ph/0302173; A.S. Jo- 

shipura and S. Mohanty, hep-ph/0302181; G. Bhattacharyya, H. Pas, L. Song, and T.J. 

Weiler, hep-ph/0302191; K. Matsuda, T. Fukuyama, and H. Nishiura, hep-ph/0302254; 

K. Cheung, hep-ph/0302265; S. Hannestad, astro-ph/0303076; O. Elgaroy and O. La- 

hav, astro-ph/0303089. 
[17] H. Fritzsch and Z.Z. Xing, Prog. Part. Nucl. Phys. 45, 1 (2000). 



FIGURES 
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FIG. 1. Illustrative dependence of mi + 1712 + 7773 on 771,3. The WMAP result sets an upper limit 
on 7773; i.e., m 3 < 0.24 eV for both 7773 > 777,2 and 7773 < 7772 cases. 
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FIG. 2. Illustrative dependence of (m) a (for a = e,fj,,r) on 771,3. Curve a: (777) e with 777,3 < m 2', 
Curve b: (777)^ « (777),- with 7773 < 7772; Curve c: (m)^ « (t77) t with 7773 > 7772; and Curve d: (m) e 
with 7773 > 777 2 . The WMAP result leads to (m) a < 0.24 eV. 
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